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Abstract 

We consider axisymmetric traveling waves propagating on the gravity-driven flow of a liquid 
down a vertical fibre. Our starting point is the two-equation model for the flow derived in the 
study by Ruyer-Quil et al. [J. Fluid Mech. 603, 431 (2008)]. The speed, amplitude and shape of 
the traveling waves are obtained for a wide range of parameters by using asymptotic analysis and 
elements from dynamical systems theory. Four different regimes are identified corresponding to the 
predominance of four different physical effects: Advection by the flow, azimuthal curvature, inertia 
and viscous dispersion. Construction of the traveling-wave branches of solutions reveals complex 
transitions from one regime to another. A phase diagram of the different regimes in the parameter 
space is oferred. 

PACS numbers: 47.15.gm, 47.20.Ma, 47.35.Fg, 47.10.Fg 
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I. INTRODUCTION 



A liquid film flowing down a vertical fibre is often encountered in a wide variety of tech- 
nological applications such as condensers, emergency cooling of nuclear fuel rods and optical 
fibre coating. It can also serve as a simple prototype for the study of wave instabilities and 
transitions in open flow hydrodynamic and other nonlinear systems. As a consequence, this 
problem has been an active topic of both experimental and theoretical research, especially 
over the past two decades. 

Experimental studies have revealed a complex wave dynamics dominated by axisymmetric 
and localised tear-drop-like structures which continuously interact with each other [1-5]. 
These structures are robust as they propagate over long distances without changing their 
speed or shape significantly; they are also separated by portions of nearly flat film, and they 
can be referred to as traveling waves. When the portions of nearly flat film between these 
structures are much longer than their characteristic length, the structures can be referred 
to as solitary waves. The formation of solitary/traveling waves results from the interplay 
between two different instability mechanisms: (i) The classical instability of a liquid film 
flowing down an inclined planar substrate prompted by inertia effects. This mode was 
initially characterised experimentally and theoretically by Kapitza and his son [6]; (ii) The 
interfacial instability of a liquid cylinder, considered first in experiments by Plateau and 
was theoretically explained by Lord Rayleigh [7, 8]. These two mechanisms are hereinafter 
referred to as the K and RP modes of instability, respectively. 

The experimental study by Duprat et al. [5], in particular, provided reports on the 
traveling waves' characteristics, namely shape, speed and amplitude, for very viscous fluids. 
The parameter values in the experiments were chosen in order to investigate the interplay 
of the K and RP modes on the waves; thus the study by Duprat et al. [5] completed 
the flow regime portrait obtained by Kliakhandler [2] for very viscous fluids but in the 
inertialess limit. Regular wavetrains were generated by means of a periodic forcing at the 
inlet. Two flow regimes were identified. For thin fibres and/or small flow rates, the RP 
mode is dominant and the solitary waves resemble beads or sliding drops whose shape is 
affected by gravity. In fact, when rescaled by the amplitude and the tail length, the profiles 
are nearly superimposed. At the same time, the flow field in a frame moving with the beads 
is characterised by recirculation zones within the beads. When closed streamlines exist in 
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the moving frame, a fluid particle is trapped in both moving and laboratory frames. Hence, 
the beads transport the trapped fluid mass downstream. For these reasons, we refer to the 
flow regime observed for thin fibres and/or small flow rates, as the drop-like regime. 

For thick fibres and/or large flow rates, the K mode is dominant and a steepening of 
the wave front is observed with an increase of the wave amplitude. Mass transport is not 
observed in this regime except for a few cases corresponding to the largest waves. We refer to 
this regime as the wave-like regime. We can conjecture that the onset of recirculation zones 
in the wave-like regime is a signature of the increased prevalence of the K mode. Eventually, 
we observe a transition from the drag-gravity regime, where inertia plays a perturbative 
role, to the drag-inertia regime, where inertia effects become dominant. Similar regimes 
and a transition between the two as inertia effects increase were first observed in the planar 
case [9, 10]. 

Noteworthy is that the drop-like regime is similar to the one observed by Quere [1] on 
a fibre or wire being pulled out form a liquid bath and results from the same physical 
mechanisms. The thin annular film that coats the wire can break up into drops. This drop 
formation process occurs when the typical time of growth of the RP instability is much 
smaller than the typical time of advection of a structure by the flow, i.e. for small fibre radii 
and/or small flow rates. 

At the theoretical front, a number of modelling approaches have been proposed within 
the framework of the long-wave approximation of the Navier-Stokes equations and associ- 
ated wall and free-surface boundary conditions [2, 3, 11-16]: The basic assumption of this 
approximation is that of slow spatial and time modulations of the film thickness motivating 
the introduction of a formal perturbation parameter, the 'long-wave' or 'film parameter' e 
measuring such modulations. Perturbation expansions in terms of this parameter then lead 
to substantial simplifications of the governing equations and boundary conditions. Addi- 
tional assumptions lead to further simplifications, i.e. small film thickness h in comparison 
to the fibre radius R or negligible inertia. The resulting models are either single evolution 
equations for the film thickness h, e.g. the model by Frenkel [11] based on the long- wave ap- 
proximation only, or systems of two coupled evolution equations for the film thickness h and 
streamwise flow rate q which combine the long-wave approximation and other approaches, 
e.g. the model by Trifonov [12] based on the 'integral-boundary-layer' approximation and 
the more recent model by Novbari and Oron [16] based on an 'energy integral' method. 
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It should be noted that all the above studies neglected the second-order viscous terms 
originating from the stream-wise momentum equation (streamwise viscous diffusion) and 
tangential stress balance (second-order contributions to the tangential stress at the free sur- 
face). The recent study by Ruyer-Quil et al. [17] formulated a two-evolution equation 
model for h and q that took into account the second-order viscous terms but also included 
inertia and was not limited to small aspect ratios h/R. The model was based on a combi- 
nation of the long-wave approximation and a weighted-residuals approach using appropriate 
polynomial test functions for the velocity field - a 'weighted residuals integral boundary 
layer' (WRIBL) model following the terminology introduced by Oron [18]. It should be 
noted that the WRIBL model is consistent [19] at 0(e) for the inertial terms and at O(e^) 
for the remaining contributions, whereas the models obtained by Trifonov [12] and Novbari 
and Oron [16] are not consistent at 0(e). Furthermore, the study by Duprat et al. [20] 
compared the wavetrains generated experimentally by periodic forcing at the inlet to the 
traveling-wave solutions of the WRIBL model showing remarkable agreement in all cases, 
thus validating experimentally the model. Experimental validation was done in the study 
by Ruyer-Quil et al. where the traveling-wave solutions of the WRIBL model were favorably 
compared to the experiments by Kliakhandler et al. [2] while the spatio-temporal dynamics 
of the film computed numerically with the WRIBL model was shown in be in good agreement 
with the experiments by Duprat et al. [4] 

As was shown in [17], the second-order viscous terms have a dispersive effect on the 
speed of the linear waves (they introduce a wavenumber dependence on the speed) and 
hence we shall refer to this effect as viscous dispersion. Viscous dispersion influences the 
shape of the capillary ripples in front of a solitary hump, more specifically, the amplitude 
and frequency of the capillary ripples, an effect which is amplified as the Reynolds number 
is increased. Hence, viscous dispersion is in fact a linear effect, but interestingly it can have 
some crucial consequences on the nonlinear dynamics of the film and the wave- selection 
process in the spatio-temporal evolution. After all, solitary pulses interact through their 
tails which overlap, i.e. the capillary ripples and their amplitude and frequency will affect 
the separation distance between the pulses. For example, smaller-amplitude ripples will allow 
for more overlap between the tails of neighboring pulses, thus decreasing their separation 
distance. These points have been analyzed in detail in the recent work by Pradas et al. [21] 
on coherent structures interaction on falling films on planar substrates and the influence of 
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viscous dispersion on interaction. Their analysis was based on the weighted-residuals models 
obtained in [22-24]. 

The main aim of the present study is to characterize theoretically the solitary/traveling 
waves propagating down the fibre within the framework of the WRIBL model. A first stab to 
the investigation of the traveling wave solutions of the WRIBL model was the recent study 
by Ruyer-Quil et al. [17, Sect. 6]. In this study traveling wave solutions were constructed 
numerically and were favorably compared to the experiments by Kliakhandler et al. [2], as 
noted earlier. Here we undertake an asymptotic analysis of the governing equations for 
solitary/traveling waves in various limiting cases, i.e. in the limits of small/large values of 
the pertinent parameters. We also obtain both numerically and asymptotically static drops 
in the drop-like regime. Furthermore, by using elements from dynamical systems theory, we 
provide a detailed and systematic parametric study of their speed, shape and amplitude, i.e. 
we construct bifurcation diagrams for their speed as a function of pertinent parameters, as 
well as obtain ranges in the parameter space for which the K or RP modes of instability, 
prompted by inertia and azimuthal curvature, respectively, are dominant. We scrutinize the 
four different regimes described earlier (drop-like, wave-like, drag-gravity, drag- inertia), and 
provide detailed phase diagrams and corresponding regime maps for very viscous and less 
viscous fluids, thus providing a deeper understanding of the problem as well as new insights 
and also completing the flow regime diagram provided in [5] for viscous fluids. 

Section II introduces the pertinent non-dimensionalization and the WRIBL model. Soli- 
tary wave solutions are constructed in § III. Asymptotic limits for small/large values of the 
pertinent parameters are analyzed in § IV. Traveling waves corresponding to the experi- 
mental conditions considered in [2, 4] are next discussed in § V. Section VI presents a phase 
diagram of the different regimes. Finally, a summary of our flndings and concluding remarks 
are offered in § VII. 

II. FORMULATION 

A. Natural set of parameters 

Consider a fllm flowing down a vertical cylinder under the action of gravity. The liquid 
has dynamic and kinematic viscosity, fj, and u, respectively, density p and surface tension 
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cr. The flow is assumed to remain axisymmetric. f, x, u and t denote the radial, pointing 
outwards from the fibre centreline coordinate, the axial coordinate oriented along gravity, the 
axial velocity distribution and time, respectively (bars are used to distinguish dimensional 
from dimensionless quantities unless the distinction is unnecessary). From simple physical 
considerations and without prior knowledge of the specific details of the system, the following 
scales can be readily identified: The fibre radius R, the Nusselt thickness /in of the uniformly 
coated film, the length and time scales, li, = v'^l'^g~^l'^ and = v^^^g'^^^, based on gravity 
and viscosity (making explicit the balance between gravity and viscous forces giving rise to 
the Nusselt fiat-film solution) and the capillary length Ic = \foJ{j)g). 

A first set of pertinent dimensionless groups arises from these scales. The aspect ratio 

a = h^/R (la) 

which assesses azimuthal curvature effects at the scale of the film, the Goucher number [25], 

Go = R/lc, (lb) 

that compares azimuthal and axial surface tension effects and the Kapitza number, 

r ^ a/{pu'/'g'/') = {IJl^f , (Ic) 

comparing surface tension and viscosity. Useful combinations of these parameters are = 
h^/lu and h^/lc- The former compares the film thickness to the gravity- viscous length scale 
and, indirectly, inertia and viscosity since the Nusselt base flow is the result of the balance 
of gravity and viscosity. The latter, h^/lc is related to the Bond number Bo = pgh^/a = 
{h^/lcY' comparing surface tension and gravity at the scale of the film. 

The advantage of the set of parameters a. Go and F is that when the geometry and 
the working fluid are fixed, the Goucher and the Kapitza numbers Go and F are constant 
and the only free parameter is a. From an experimental point of view, a can be varied 
independently by varying the inlet flow rate. The Kapitza number F is entirely defined by 
the fluid properties independently of the flow characteristics, whereas the Goucher number 
Go can be easily varied by replacing the fibre. Hence, the parameters a. Go and F can 
therefore be viewed as 'natural' for the fibre problem, and we will systematically recast our 
results in terms of these parameters a. Go and F. 

Table I gives the physical properties of four different fluids of increasing viscosities com- 
monly used in experiments and corresponding to a wide range of Kapitza numbers together 
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V (mm^s 


p (kgm ^) 


o (mNm 


If. (mm) / 


y (mm) 


r 


water 


1 


998 


72.5 


2.7 


0.047 


3376 


Rhodorsil silicon oil v50 


50 


963 


20.8 


1.5 


0.63 


5.48 


castor oil 


440 


961 


31 


1.8 


2.7 


0.45 


Rhodorsil silicon oil vlOOO 


1000 


980 


21.1 


1.5 


4.7 


0.10 



TABLE I. Fluid properties, capillary length 1^^ gravity-viscous length ly and Kapitza number used 
in the present study. The data for silicon oil v50 and castor oil have been taken from [2, 4]. 

with the corresponding capillary lengths and viscous-gravity length scales. For simplicity, 
our results will be presented for the Kapitza numbers listed in table I. 



B. WRIBL model 

We now adapt Shkadov's scaling [26] and introduce different length scales for the stream- 
wise (axial) and radial directions. The length scale in the radial direction is the Nusselt 
thickness /in, whereas the length scale in the streamwwise direction is chosen as k/in de- 
fined by the balance of the streamwise pressure gradient induced by capillarity, oc (jd^xxh, 
and gravity acceleration, pg, which gives k = [a / {pgh'^)Y^^ = {Ic/h^Y/'^. The time scale 
is defined with reference to the Nusselt solution of uniform thickness (a result of the bal- 
ance of gravity and viscosity). The volumetric fiow rate per unit length of circumference, 
gN = ju^^^ urdr, of a film of constant thickness /in is given by 

g-N = ^0(«) , (2) 

where is a geometric factor defined in (Ala) and measures the deviation of the fiow- 
rate-to-thickness relation from the cubic dependency corresponding to the planar geometry 
(0(0) = 1). Similarly to the streamwise length scale, the time scale is stretched by a factor 
K and thus defined as ?>tdi^/q^ = un/lght^cp^a)]. 

Our basic equations for the analysis to follow are the WRIBL model obtained in [17], 
a set of two evolution equations for the local film thickness h{x,t) and the local fiow rate 
q[x, t) = jR+hi^,t) ^ ^^^^ Pqj, ^j^g gg^j^g q£ clarity and completeness we re- write the WRIBL 
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model here, 
dth 
6dtq 



ah 



(3a) 



-F{ah) j^dr^q + G{ah) j^dxh 
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+7] 



I{ah) 



(a 



30(a) q 
(f){ah) 



J{ah) — {dxhY — K{ah)-^—^ L{ah)—dxxh + M{ah)dxxq 



(3b) 



in terms of the Shkadov scahng, where 0, F, G, /, J, L, and M are functions of 
the aspect ratio a given in Appendix A. Equation (3a) is the (exact) dimensionless mass 
balance whereas (3b) is the streamwise momentum equation averaged across the film with 
a weighted-residuals approach. It should be emphasized that the WRIBL model has been 
validated in [5, 17] through direct comparisons to the experiments in [2, 4, 5] as noted in 
§ I [for both very viscous and less viscous liquids (castor oil and silicon oil V50, see table I) 
and a wide range of the parameters (0.15 < Go < 1 and 0.5 < a < 4.5)]. 

Shkadov's scales introduce three new dimensionless groups besides the aspect ratio a = 
h]<f/R, a 'reduced Reynolds number', 



S = Sq^^s/i^un) = (aGo) 



11/3 



(a 



)r3/2, 



(4a) 



which compares inertia and the viscous drag at the scale k/in introduced by the balance of 
gravity and capillarity, a streamwise 'viscous dispersion parameter'. 



\4/3 



and the dimensionless group. 



(3 = a'/f] = ~o?l''Go-'l\ 



(4b) 



(4c) 



which is a combination of a and ?7 and compares azimuthal to axial surface tension effects. 
We have made explicit in (4) the relations of 5, 77 and /3 to the 'natural' parameters a. Go 
and r. Notice that all second-order /viscous-dispersion terms are gathered in the last row of 
(3b) and are multiplied by 77. Finally, we recall the usual definition of the Reynolds number 
based on the flow rate. It = q^jv = h^(j){a)/?> where again /i^ — h'^^/ly. 

The advantage of Shkadov's scaling stems from (i) the direct reference to the Nusselt 
uniform film flow that simplifies the comparisons between solutions, with the Nusselt solution 
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corresponding to constant values of the film thickness and flow rates h = 1 and q = 1/3; 
(ii) the association of a single parameter to each physical effect affecting the balance of 
the different forces: Inertia {6), azimuthal surface tension viscous dispersion (77) and 
geometry (a). 

As noted in § I the WRIBL model was derived with the long-wave approximation (i.e. 
under the assumption of slow space and time modulations, d^^t ~ where e is the long- 
wave/film parameter) and a weighted-residuals approach in which the velocity field is ex- 
panded on an appropriately chosen set of test functions. This expansion takes into account 
the (small) deviations of the velocity field from the uniform-thickness solution. As also 
noted in § I, contrary to the model obtained by Trifonov [12] - see also [27] - and to the 
model by Novbari and Oron [16], (3) is consistent up to 0(e) for the inertia terms and up to 
O(e^) for the remaining contributions (and accounts for viscous dispersion). Indeed, both 
Trifonov's and Novbari and Oron's approaches assume a self-similar velocity distribution 
and do not account for the deviations of the velocity profile induced by the free-surface 
deformations. For this reason, their two-equation formulations lack consistency even at first 
order in the film parameter. Furthermore, we note that the energy-integral approach em- 
ployed by Novbari and Oron [16] is not consistent with the kinetic energy balance of the 
flow. Indeed, writing formally the axial momentum equation as Ai{u) = 0, Novbari and 
Oron's averaged momentum equation reads j^^^ M.{u)udr = whereas the kinetic energy 
balance of a section of the liquid corresponds to j^^^ M.{u)ud{r'^) = 0. Truncating then 
Ai{u) at O(e^) is equivalent to the Galerkin approach that can be used to reduce the algebra 
leading to (3b) [17]. Noteworthy is that the two-equation model (3) is not limited to small 
aspect ratios unlike, e.g. the model by Roberts and Li [15]. 

To end this section let us point out one apparent drawback of the Shkadov scaling, namely 
the divergence of the dimensionless parameter /3 = a'^/t] as the viscous dispersion parameter 
goes to zero for a = 0(1). The divergence of /3 signals that the typical length of a wave 
is not determined by the balance of gravitational and axial capillary forces, as assumed in 
Shkadov's scaling, but rather by the balance of axial and azimuthal capillary forces, in which 
case the typical curvatures of the beads in the azimuthal and axial directions must coincide. 
The beads have thus a drop-like rounded shape, the long-wave approximation starts to be 
violated and the viscous dispersion effects cannot be a priori discarded as in the planar case. 
Considering drop- like beads (see §-IVA) it can thus be useful to adopt a scaling based on 
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the radius R of the fibre, which gives the time scale v/ (gR). This scahng introduces a Gahlei 
number, Ga = gR? /v"^ = Gb^F^/^. 



C. Surface equations and saturation numbers 

Neglecting inertia and viscous dispersion, Craster and Matar [3] formulated a single 
evolution equation for the film thickness h, the Craster and Matar (CM) equation. 



h'<l>{&h) f 13 



2 

For sufficiently thin films, that is a — !■ 0, we obtain from (5) 



0. (5) 



dth + 



0, (6) 



which is the equation derived initially by Frenkel [11] (see also [13]). Equations (5) and (6) 
are the simplest evolution equations balancing all relevant physical effects, gravity, viscous 
drag, surface tension and the fibre curvature. They are equations for the free surface h{x, t) 
only and following the terminology introduced by Ooshida [9], we shall refer to them as 
surface equations. 

The CM equation offers a simple prototype for the flow, easily amenable to mathematical 
and numerical scrutiny. It can be obtained asymptotically from (3) in the limit 5 — ?■ 
and ?7 — !■ when the nonlinear term —l/2dx (<5/(l + ah){dxKf') of the azimuthal curvature 
gradient is omitted. Yet, as already stated, a = 0(1) implies that /3 diverges to inflnity and 
therefore that second-order viscous dispersive terms cannot be a priori discarded and that 
the long-wave assumption is invalid. By contrast, the derivation of the Frenkel's equation 
(6) can be obtained in the distinguished limit a — ?■ 0, r/ — and /3 = a^/rj = 0(1). 
Nevertheless, comparisons of the solutions to the CM equation to the experiments show 
good agreement [17, 20]. 

Using the surface equation (6), Kalliadasis and Chang [13], and Chang and Demekhin [28] 
analysed the mechanism of the formation of drops observed by Qyere [1] when a flbre or wire 
is drawn out of a liquid bath. Quere observed suppression of the RP mode for sufficiently 
small coating films. More specifically, below a certain critical thickness the film deposited 
on the wire develops small-amplitude interfacial waves with the fiow preventing their growth 
into drops (such drops would always be observed when the wire is horizontal, as the sup- 
pression of the RP mode induced by the fiow is absent in this case). On the other hand, 
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beyond the critical film thickness growth of the interfacial waves into drops was observed. 
Kalliadasis and Chang [13] found that the amplitude of the solitary-wave solutions to (6) 
diverges/'blows-up' for /3 larger than a critical value (3c = 1.413, which closely corresponds to 
the formation of drops in Quere's experiments. These authors also performed an analytical 
construction of the solitary waves for fi ^ using matched asymptotic expansions. They 
showed that drop formation results from the inability of the flow advection to saturate the 
instability growth through a nonlinear saturation mechanism. 

The advection and the growth of the instability can be compared through the definition 
of the advection time of an interfacial structure over its length and the definition of a 
typical time of growth of the structure Tg as the inverse of the maximum growth rate. Based 
on the Fenkel equation (6), the stability analysis of the uniform film leads to the dispersion 
relation 

c. = fc + ^(/3-P) (7) 

which governs infinitesimal perturbations around the base Nusselt flow of wavenumber k 
and angular frequency uo. Based on (7) the ratio Ta/rg reads: 

ra/Tg = UJi/uJrlr^^(xl3'"^ (8) 



(see [17] for details). Therefore, (5 compares Tq to Tg. For, /3 < /3c, the instability growth is 
slower and thwarted than the advection by the flow. The same mechanism is also in play in 
the saturation of the drops though it is then strongly nonlinear. For these reasons, we refer 
to /3 as a saturation number., a term that was first introduced in [5]. 

From the CM equation (5), we get Ta/Tg oc {(3*)'^^'^ where the composite parameter (3* is 
defined as [4, 17]: 

r = /3c-'/^(l + «)-«/^ (9) 

Cfc is the speed at which infinitesimal deformations of the free surface are kinematically 
transported by the flow, or kinematic wave speed: 

1 



Ck 



a(f)'(a) 
30(a) 



8(fe-l)(21og(&)6^-&^ + l) 
3(41og(6)64-364 + 462_i) ' ^ 



. + a 
with b = 1 + a. 

Since Ck{a = 0) = 1, limQ,_j.o /?* = (3. In this limit the Frenkel equation (6) limited to 
very thin films (a ^ 1) applies and /3* is a generalisation of the saturation number /3 to film 
flows with thicknesses comparable to the fibre radius. 
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III. SOLITARY WAVES AND DYNAMICAL SYSTEMS THEORY 



Experimental studies [2-5] reported the formation of axisymmetric traveling waves (TWs) 
propagating without deformations and at constant speed over long distances. Solitary waves, 
i.e. TWs separated by constant-thickness layers of fluid, or substrates, much longer than the 
characteristic length of the waves, were commonly observed sufficiently far downstream. 
Theoretically, solitary waves can be viewed as periodic TWs with an infinitely long wave- 
length. The aim of this section is to investigate infinite-domain solitary waves using elements 
from dynamical systems theory. 

In a frame of reference moving with the speed c of the waves, ^ = x — ct, the fiow is 
stationary and the set of partial differential equations (3) can be converted into a set of 
ordinary differential ones. The mass balance (3a) can be integrated once. 



a 



h+ -h' 



(11) 



where go is the rate at which the fluid flows under the waves. The averaged momentum 
balance (3b) next reads. 
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cq' -Fi&h) ^q' + G(ah)^h' 
h h'' 



+ 



I{dih) 



[a 



30(a) q 
(f){ah) 



h{l + h'" + 



/3 



-7] 



[l + ah) 



a 



J{ah)f^{h'f 



K{ah) 



1 + ah 
q'h' 



h' 



h" - 



a 



l\2 



21 + ah 



h 



Lidih)-h" + Miah)q" 
h 



0, 



(12) 



where the primes denote differentiation with respect to the moving coordinate ^. Using (11), 
(12) can be recast into h'" = f{h, h', h") , where / is a function of the thickness /i, its flrst 
and second derivatives and the parameters 5, a, rj and c. We thus end up with a dynamical 
system of dimension three, 

d 
di 

where U = (f/i, U2, U^Y = (/i, h', h'J. 



-V = iU2,UsJiUi,U2,U3)y 



(13) 



Solitary waves correspond to homoclinic orbits in the phase space connecting a flxed 
point to itself. Here we restrict our attention to single-loop homoclinic orbits corresponding 
to single-hump solitary waves in real space. The flxed points of the dynamical system (13) 
satisfy h' = h" = and, 

h^(j){ah) f a 2\ f..s 
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Li 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

C/Cfc 

FIG. 1. Film thicknesses hi = 1 and hu corresponding to the location of the fixed points as function 
of the ratio of the normalised wave speed c to the kinematic wave speed defined in (10). Solid, 
dashed and dotted lines refer to a = (planar case), a = 0.5 and a = 1, respectively. 

Requiring that /i = 1 is a solution to (14) gives: 

+ (15) 

In addition to the solution h = 1, there is one more real positive solution to (14) with (15) 
and the dynamical system (13) admits two fixed points Ui = (1,0,0)* and Un = (/in, 0,0)* 
whose positions are displayed in Fig. 1 as functions of the wave speed c and aspect ratio a. 
The two fixed points coincide when the wave speed c is equal to the speed of linear waves 
of infinitesimal amplitude and infinite length that are neutrally stable. This situation corre- 
sponds to the definition of the linear kinematic waves in the zero-wavenumber limit, whose 
speed is given in (10). In the limit a = 0, we recover the expression hu = — l/2+-\/3(c— 1/4) 
corresponding to a film fiowing down an inclined plane [10, 29]. We note that extending 
this expression to the axisymmetric geometry by replacing c with c/c^, gives a reasonable 
approximation hu ~ — 1/2 + -\/3[(c/cfc) — 1/4] for the position of the second fixed point even 
for a = 0(1) (see Fig. 1). Finally, we note that it is sufficient to consider homoclinic orbits 
around only one of the two fixed points because of the presence of a continuous family of 
Nusselt fiat-film solutions parameterized by the reduced Reynolds number 6 (or h-^ = h^/ly) 
when Go and F are held constant. Indeed, homoclinic orbits connecting Un correspond to 
phase-space trajectories connecting Ui through the transformation h^ h^hu. 

The shape of the tail and front of a solitary wave can be determined by considering how 
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the corresponding homoclinic orbit in the phase space approaches and leaves the fixed point 
it connects to. Let us consider the hnear stabihty of the fixed point Ui. The dispersion 
relation governing infinitesimal perturbations ~ exp(A^) is 



+ X^riDr, + XSDs - 3(1 + a) 



Ck) 



0, 



(16) 



where 



a)M 



and Df, 



G_ 
9" 



+ 



/3 



5(l + a)2 



Equation (16) can be reduced to the canonical form P{y) 
of variable X = y — t]D^/3, where 
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+ py + q = hj the change 
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Using the Cardan formulae, (16) admits a real eigenvalue and a complex conjugate pair 
when the discriminant A = Ap^ + 21q^ is > 0. When A < 0, (16) admits three real 
eigenvalues. Therefore Ui changes from a saddle to a saddle-focus at A = 0. When Ui is 
a saddle, the homoclinic orbit departs and returns to the fixed point monotonically along 
the two eigenspaces corresponding to the eigenvalues of smallest absolute value, and the 
corresponding tail and front of the solitary wave are monotonic. Conversely, when Ui is a 
saddle-focus, the homoclinic orbit leaves monotonically along the eigenspace spanned by the 
eigenvector corresponding to the real eigenvalue and returns to the fixed point by spiralling 
on the eigenspace spanned by the eigenvectors corresponding to the complex pair. This 
spiral corresponds to capillary ripples at the front of the solitary wave while the tail remains 
monotonic. 

Figure 2 displays the speed and maximum amplitude of the solitary waves corresponding 
to homoclinic orbits around Uj = (1,0,0)*. Parameters correspond to the fiuid properties 
of silicon oil v50 (F = 5.48, see table I) and to different fibre radii, R = 1.5 mm, 0.35 mm 
and 0.25 mm [Go = 1, 0.24 and 0.17). The solutions have been computed by continuation 
using AuT097 together with its subroutine HomCont [30]. In all cases A > 0, i.e. Ui is 
a saddle-focus and the solitary waves are characterised by monotonic tails and oscillatory 
fronts. 

Two different behaviours can be observed for small and large thicknesses, or for S <^ 1 
and (5^1. [31] For small thicknesses {6 ^ 1), the maximum height /imax and speed of the 
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FIG. 2. Speed c (left) and maximum height /imax (right) of the solitary waves as fmiction of the 
reduced Reynolds number 6 for different fibre radii: R = 1.5 mm (Curves 1), 0.35 mm (2) and 
0.25 mm (3). Solutions to (3) (solid lines) are compared to the solutions of the CM equation (5) 
and of the model (3) when inertial terms are set to zero (dashed and dashed-dotted lines). The 
fluid properties correspond to Rhodorsil silicon oil v50 (see table 1). 

solitary waves exhibit local maxima which strongly increase for thin fibres corresponding 
to a stronger RP instability. In Figs. 2(a) and 2(b), the characteristics of the solitary 
wave solutions of the WRIBL model (3) and of the CM equation (5) are compared showing 
reasonable agreement. As the CM equation is parameterized by the aspect ratio a and /3, 
the local maxima are the result of the balance of curvature effects and the advection by 
the flow. This is reminiscent of the sharp increase, or 'blow-up', of speed and amplitude 
observed by Kalliadasis and Chang [13] at /3 = /3c ~ 1.413 in their study of the solitary wave 
solutions of the Frenkel equation (6). 

One might expect that the sharp increase of the local maxima of the speed and amplitude 
observed by lowering the fibre radius is related to the nonlinear saturation mechanism of 
the instability by the advection discussed earlier and, therefore, it should be correlated with 
the saturation number The validity of this hypothesis is checked in Fig. 3 where the 
maximum height of the solitary waves is depicted as a function of /3*. At a given value of 
the Goucher number Go, /3* reaches a maximum for a ~ 0.44 and tends to zero for a — )■ 
and a — > oo. For a fixed radius R, 6 and a have the same trend. This explains why /3* tends 
to zero when 6^0 and 5 — )■ oo and justifies the shapes of the curves. The local maximum 
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FIG. 3. Maximum height of sohtary waves versus saturation parameter 13* . See also caption of 
Fig. 2. 

of /imax occurs at /3* close to the maximum reached by this parameter as 5 is varied. The 
increase of the local maximum of /imax is related to the increase of the maximal value of (3* 
achieved as the ratio Go is lowered. 

We have also computed solutions of model (3) when the inertial terms are cancelled 
(5 — 7- 0). The results are shown as dashed-dotted lines in Figs. 2(a) and 2(b). As they 
asymptote to the solid lines corresponding to (3) in the limit 5 ^ 1, we can conclude that 
the difference in speed and amplitude between solutions to (3) and to the CM equation 
result from the viscous dispersion effects that contribute to the increase of the speed and 
amplitude of the solitary waves. 

Figure 2 reveals a sharp increase of the maximum height and speed of the waves above 
(5^2 corresponding to the predominance of the K mode (as the planar case is approached, 
i.e. when the fibre radius increases, this sharp increase corresponds precisely to the transition 
between the drag-gravity and drag-inertia regimes). The characteristics of the waves in this 
region will be considered in detail later on in §IVB. The separation between the local 
maxima for the speed and amplitude corresponding to the RP-dominated waves (5^1) 
and the K-dominated waves at large 6 increases for even more viscous fluids like the castor 
oil used by Kliakhandler [2] (cf. Table I) as can be observed from Fig. 4(b) where the 
solutions of the model (3) with and without inertia, and of the CM equation are compared 
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FIG. 4. Speed c (left) and maximum height /imax (right) as function of the reduced Reynolds 
number 6 for different fibre radii: R = 1.83 mm (Curves 1), 0.43 mm (2) and 0.31 mm (3). 
Solutions to (3) (solid lines) are compared to the solutions of the CM equation (5) and of the 
model (3) when inertial terms are set to zero (dashed and dashed-dotted lines). Homoclinic orbits 
connecting a saddle-focus (saddle) fixed point to itself are shown in thick (thin) solid lines. Squares 
indicate the loci of saddle to saddle-focus transitions. The fluid properties correspond to castor oil 
(r = 0.45). 

for fibre radii corresponding to the same ratios Go = 1, 0.24 and 0.17 as in figure 2. This 
increase can be understood by considering the definition of /?*, which is a function of the 
aspect ratio a and Go, and the definition of 5 = (aGb)^^/^0(a)r^/^. Since the Kapitza 
number T decreases with the kinematic viscosity, the maximum of /3* for a given value of 
the Goucher number Go corresponds to smaller and smaller values of the reduced Reynolds 
number 6 as the viscosity of the fluid is increased. In contrast with the results for silicon oil 
v50 (cf. Fig. 2), a transition from a saddle-focus to a saddle fixed point has been observed 
corresponding to the disappearance of the capillary ripples at the front of the solitary wave 
at values of 6 above 0.1. The precise location of this transition varies with the fibre radius 
R and is indicated by squares in figure 4(a). 

Conversely, the separation of the solitary wave characteristics, such as speed and maxi- 
mum height, as a function of 6 into two distinct regions, at low and high reduced Reynolds 
numbers, vanishes at low viscosities. Indeed at low viscosities, or equivalently, high Kapitza 
numbers, the RP mode occurs at relatively high values of 6 where the K mode already takes 
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FIG. 5. (a) Maximum height /imax of the sohtary waves as function of 6. Solutions of (3) (sohd 
hues) are compared to the solutions of the CM equation (5) (dashed lines) and solutions of (3) 
where the inertia terms are set to zero (dashed-dotted lines), (b) Maximum height /imax (solid 
line) and parameter /3* (dashed lines). Parameters correspond to water (F = 3376) and different 
fibre radii: R = 2.75 mm (Curves 1), 0.64 mm (2) and 0.46 mm (3). 



over. In Fig. 5(a) we have redrawn Fig. 2(b) for water (F = 3376, see Table I), which is 
fifty times less viscous than silicon oil v50. The figure compares the maximum height of the 
solitary waves for a ratio of the Goucher number Go equal to the ones chosen for the compu- 
tations shown in Figs. 2, to the amplitude of the solutions of the CM equation (5). At small 
values of 6, the amplitude of the waves is significantly larger than the amplitude of the solu- 
tions of the CM equation, which signals the influence of the K mode on the RP instability. 
Notice that this effect cannot arise from the second-order viscous terms as the cancellation 
of the inertial terms leads to results comparable to the solutions of the CM equation. On 
the other hand, the influence of the RP instability on the K mode is illustrated in Fig. 5(b) 
where /3* and /imax are plotted versus 6. At small fibre radii R = 0.64 mm and R = 0.46 mm 
{Go = 0.24 and 0.17), the characteristic sharp increase of the solitary-wave amplitude is dis- 
placed to values of 6 smaller than 6^2 corresponding to a generalised saturation number 
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^111 ax 
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Cu/{gR^) Cdrops 1^/(5^^) 


0.236 


0.23 
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3.1 


0.77 


0.22 


0.25 


0.168 


0.15 
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0.30 


0.110 


0.075 


18 


38.5 


1.36 


0.24 


0.29 


0.055 


0.022 


81.2 


383 


1.80 


0.19 


0.22 


0.044 


0.014 


129 


773 


1.85 


0.16 


0.18 



TABLE II. Local maxima of the speed c and amplitude /imax with respect to 5 for given values of 
the Goucher number Go, or equivalently a, obtained from the CM equation (5). Cdrops refers to the 
speed of quasi-static drops sliding downwards coated fibres (see text and Appendix B). 

IV. LIMITING CASES 

In this section, we focus on the different regions of the wave-characteristics' diagrams 
displayed in Figs. 2, 4 and 5, and consider all possible asymptotic limits. 

A. Small Goucher number limit: the drop-like regime 

Let us first consider the local maxima of the amplitude /imax and speed c with respect to 
5 for given values of the Goucher number Gb, or equivalently a, observed for viscous fluids 
like silicon oils in the inertia-less limit (5 ^ 1). Table II depicts these quantities as obtained 
from the CM equation. As Go tends to zero, the RP instability mechanism becomes more 
and more efficient and we observe a sharp increase of the local maxima of the amplitude /imax 
and of the maxima of the speed of the waves. For such waves, the amplitude can be several 
orders of magnitude larger than the Nusselt fiat film on which they stand, and variations of 
the Nusselt thickness should only slightly modify the wave characteristics (except perhaps 
when the film becomes so thin that the corner dissipation at the wave front dominates over 
the dissipation in the bulk). 

Since R <^ /c, azimuthal surface tension effects dominate over gravity and the typical 
length scale and height of a wave should correspond to the radius R of the fibre as already 
pointed out in § II B. The wave speed should then be determined by the balance of viscosity 
and gravity at the scale R which gives a typical velocity of gR? /v for viscous-gravitational 
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drainage. Justification of tfie neglect of the inertial terms demands that the Gahlei number 
Ch = Gb'^r^/^ is small, which is satisfied for all tested fluids except for water (F = 3376). 
Our computations of the solutions to the CM equation (5) confirm these scaling arguments 
as shown in table II. 

In Fig. 6, we contrast the wave profiles rescaled with respect to R corresponding to 
Table II. Except from the front and back of the waves corresponding to the return to the 
fixed point, the wave profiles are rather symmetric. This front-to-back symmetry shows that 
gravity does not affect the wave profile, as expected, since the typical size of the wave R is 
much smaller than the capillary length 1^. Therefore, solitary waves in this regime resemble 
isolated drops sliding under the action of gravity on a wettable fibre, which is precisely why 
we refer to this regime as the drop-like regime. [When R is larger and approaches Ic, on the 
other hand, the drops will feel the effect of gravity as they grow and they will eventually 
resemble falling pendant drops.] It corresponds to the observation by Quere [1, 25] in the 
coating of wires or thin fibres drawn out of a bath of viscous liquids that the thin annular 
film deposited on the wires/fibres breaks up into drops. 

We have checked this analogy by computing the shape of static drops with zero con- 
tact angles sitting on a fibre coated with a substrate film of the same liquid (details of the 
calculation are given in Appendix B). The agreement of the wave shapes to the symmet- 
rical static-drop shapes is remarkable even in the case of nearly spherical drops such as 
those obtained at Gb = 0.044, where the long-wave assumption is no longer valid. We can 
therefore conclude that the CM equation (5), and by extension the WRIBL model (3), are 
accurate in the drop-like regime where surface tension is dominant and even if the long- wave 
approximation does not hold. Besides, the remarkable agreement between the solutions to 
the CM equation (5) and the WRIBL model (3) already noted in Fig. 2, shows that the 
second-order streamwise viscous terms do not affect the amplitude and speed of the drops. 

Following [13], an analytical estimate of the amplitude and speed of the drop-like waves 
in the limit Go ^ may be obtained via matched asymptotic expansions. The appropriate 
small parameter is the dimensionless speed of the drops. By balancing viscous and capillary 
forces at the back of the waves, one can easily extend to sliding drops the Landau-Levich- 
Derjaguin law obtained by Quere [25] in the case of fibres drawn out of a bath. The speed of 
sliding drops is thus governed by Eq. (B4) which compares favorably to the results from the 
CM equation in table II. As a matter of fact, this agreement shows that the thickness of the 
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FIG. 6. Wave profiles corresponding to the solutions to the CM equation (solid lines) and static 
drop shape (thin dashed lines). Labels correspond to the Goucher number. Values of the other 
parameters are listed in table II. 

residual film on which the drops slide is determined by the balance of surface tension and 
viscous dissipation in the meniscus region. To estimate the speed and amplitude of drop- 
like waves, one would have to take into account the gravity acceleration and higher-order 
corrections in the outer region, namely the viscous dissipation in the drop. A task which 
is difficult, as (a) with the Frenkel equation resolving fully the leading-order outer region 
requires matching up to third order [13] and quite likely this is the case here, (b) a single 
'composite equation' for the whole domain, i.e. for both the drop and residuals films, as e.g. 
in the 'drag-out' problem in coating theory [32], does not exist. 

B. Large 5 limit: the drag-inertia regime 

To understand the change of behaviour of the solitary- wave characteristics for 5^1, 
we look at the wave profiles and their projections on the plane [h, h'). Figure 7 compares 
two single-loop homoclinic orbits for a small and a large value of 6. For 5^1, solitary 
waves have a relatively symmetric shape. Except from the neighbourhood of the fixed point 
Ui, where the escape from Ui is monotonic and the return to it is oscillatory, a symmetry 
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FIG. 7. Profile (a) and projected trajectory (b) onto the plane (/i, h') of single-loop homoclinic 
solutions to (3) corresponding to single-hump solitary waves. Parameters correspond to Rhodorsil 
silicon oil v50 (F = 5.48) and R = 0.25 mm. Solid and dashed lines refer to 5 = 5 and 5 = 3 10~^, 
respectively. 

between the front and back of the waves is observed with a steep front and a steep back. On 
the contrary, at large values of 5, the front and back of the solitary waves present radically 
different shapes, with a gentle sloping back edge and a steep front edge. 

The observed difference between front and back of the solitary waves in the large-(5 case 
can be explained by examining the linearised dynamics around Uj. The dispersion relation 
governing infinitesimal perturbations varying as exp(A.^) is given in (16). At a given radius 
R and capillary length Z^, large thickness, /in ^ 1, corresponds to 5 ~ /ij/^^0(a) ^ 1, 
(5/5 = {Ic/ Ry / (Slh) ^ 1 and possibly large viscous dispersion since rj ~ h^^. Let Ai be 
the positive eigenvalue corresponding to the unstable manifold of Ui. The eigenvalues of 
the tangent subspace to the unstable manifold satisfy Re(A3) < Re(A2) < 0. When Ui is 
a saddle-focus (A > 0), we further denote A2,3 by — S ± if2 with S > and > 0. From 
(16) we obtain the estimate Ai ~ 6^^ ^ 1. Since Ai + A2 + A3 = —rjDr,, we immediately 
get an estimate of the mean value, (A2 + A3)/2 ~ —77, so that S ~ 77 when A > 0. As a 
consequence, at the back of a solitary wave, the monotonic escape from the fixed point is 
slow, whereas at the front, the return to Ui is fast. The above estimates have been confirmed 
by computations of Ai and S for the solitary waves shown in Fig. 2 corresponding to silicon 
oil v50. 

Focusing now at the back of the solitary wave and defining a slow variable ^ = ^/S, (12) 
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reads, 



f3I{ah) 



h 

I{ah) 



5(1 + ahy(j){a) J de 



dh 



30(a) q 
(j){ah) 



- h 



where q is given by (11) and (15). Equation (19) can be formally rewritten as, 

g{h, c- a, /3/(5)i = -HQi, c; a) + 0{7]/5'', 5-^) , 



(19) 



(20) 



expressing the balance at the back of the solitary waves of inertia (at the left-hand side), 
viscous drag and gravity acceleration (at the right-hand side). As a consequence, the roots 
of the right-hand side of (19) correspond to the fixed points of (13). 

As the homoclinic orbit departs from Ui, h increases up to hu which is larger than unity 
since c > (cf. Fig. 1). At this point, h goes through a maximum if Q is nonzero. The 
resulting orbit is thus a heteroclinic one linking Uj and Un which contradicts the fact that 
we are considering a single-loop homoclinic orbit [10]. As a consequence, Q must go to zero 
aX h = hu which signals a 'critical film thickness' he at which inertial terms must go to zero. 
In the limit 5 — ?■ oo, we thus obtain the condition 



hii{c] a) = hc{c] a,/3/S) 



(21) 



which gives a unique solution Ccrit and then hcrit as function of a and f3/S. As the limit 
speed Ccrit is governed by the balance of inertia, wall friction and gravity acceleration, we 
may refer to this situation as the drag-inertia regime [9]. 

Having shown that Q possesses at least one root, Q can be factorised into Q = {l + ah)[c — 
Cd+iah, q/h, f3/6)][c — Cd-{ah, q/h, (3/6)] where. 



Cd±[ah, q/h) = - — - — ± 



/A 



ah,q/h,l3/5 



and A 



ah,q/h,l3/5 



hJ 



F{ah)^ G{ah) 



f3I{&h) 



5(1 + a/i)3(/.(a) 



(22a) 
(22b) 



4 l + ah_ 

are the speeds of linear dynamic waves for a uniform layer of thickness h and averaged speed 
q/h [17]. In other words, the position of the second fixed point must coincide with the 
critical layer he at which the speed of the solitary wave c is equal to the speed of one of the 
dynamic waves Cd± - in fact the fastest one with speed Cd+ which separates the flow into a 
'subcritical region' (c < q+) and a 'supercritical region' (c > Cd+). The condition hu = hcnt 
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is similar to the 'Thomas condition' derived in the mathematical treatment of periodic bores 
on steep slopes, or roll-waves [33, 34] made of the regular succession of laminar flows and 
hydraulic jumps. For 5^1, the front of the solitary waves, where surface tension arrests 
the breaking of the waves, plays a role similar to that of a hydraulic jump connecting the 
subcritical and supercritical regions of a roll wave. For this reason, we may also refer to this 
regime as the roll-wave regime. 

Considering now large but finite values of S, condition (21) is not verified and one has to 
go back to (20). At the critical point h = h^vfe have Q{hc, c; a) = and a Taylor expansion 
close to criticality gives 

he - /icrit ~ - [dcQ/ dhQ] (/icrit, Ccrit) (c - Ccrit) (23) 

whenever ^^^(/icrit, Ccrit) 7^ 0. Proceeding next to a Taylor expansion of H{h,c), (20) yields 



OhQ 



where the constants and K^t are functions of /icrit? Ccrit and the derivatives h'^T-it, h'^j-it and 
at h"^.-^^ at the critical thickness of the asymptotic solution. Since for h^q ^ 1, 6 ~ hj^ and 
V ~ ^n'^ asymptotically dominant term in the right-hand-side of (24) are the viscous 
terms Kjjr]6~'^. Consequently, the convergence of the speed of the waves to the asymptotic 
value Ccrit (a) satisfies 

c - Ccrit(«, P/S) oc r]/6^ ~ (25) 

in agreement with our numerical computations. 

Figure 8 compares the phase speed c of the solitary waves to the asymptotic limit Ccrit 
as function of 6 for the four different fluids of increasing viscosity that we considered (fluid 
properties are gathered in Table II). For weakly viscous flows, the RP and K mode reinforce 
each other and the speed c of the waves is much higher than the asymptote Ccrit (see Fig. 8(a)). 
At larger viscosities (Figs. 8(b) and 8(c)) the K mode ceases to be affected by the RP 
instability and a rapid convergence to Ccrit is observed as 6 is increased. Whereas for very 
viscous fluids like silicon oil vlOOO (one thousand times more viscous than water, cf. Table I), 
the curves start again to separate signalling the arrest of the convergence to the asymptote 
Ccrit by the axial viscous effects. 
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FIG. 8. Speed c of the solitary waves as function of the reduced Reynolds number S for three 
values of the Goucher number Go = 1.01 (curves 1), 0.236 (curves 2) and 0.168 (curves 3) for 
fluids of decreasing Kapitza numbers. Solutions to (3) (solid lines) are compared to the asymptotic 
predictions Ccrit (dashed lines). 



C. Weakly nonlinear analysis 

Let us now consider the regions of the speed and maximum height diagrams in Figs. 2, 
4 and 5 corresponding to the transition between the drop-like regime and the drag-inertia 
regime. They correspond to situations where neither the K nor the RP instability mecha- 
nisms are strong enough to promote large amplitude waves, i.e. when 5 ~ 1 and /3* ^ 1. 
We can then proceed to a weakly nonlinear analysis to characterize the shape, speed and 
amplitude of the waves. Since 5 ^ 1, we substitute h = 1 + H and q = 1/3 + Q where 
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H = 0{e) and Q = 0{e) are small deviations from the base state {h,q) = (1,1/3), (i.e. 
e ^ 1). Taking the distinguished limit e <^ e^, we are led from the WRIBL (3) model to a 
single evolution equation for the deviation H: 



dtH + c,d,H + ^a)Hd,H + 

' (a 



/3 



+6 



3J(a) 



+r] 



dttH 



[a) 



3(l + a)3 
F(a 



dxxH + 
^-dtxH^ 



1 



3(1 + a 
G{a) 



9(1 + a 
Lia) 



dxxxxH 
-dxxH 



3(1 + a) 



where 



$(a) 



3(2 + 5)0 + a[(6 + 5a)0' + a(l + a)(p"] 



0(e6^eV), (26) 



(27) 



3(1 +a)20 

is a function of the aspect ratio a. Noteworthy is that a stability analysis of the base state 
(/i, q) = (1, 1/3) based on (26) leads to the same dispersion relation as for (3). 



1. Drag- gravity regime 

We consider here the limit of a thin film compared to the fibre radius, a <^ 1, when 
the radius of the fibre is constant or equivalently Go is constant. Since r] = {h-n/lcY^^ = 
viscous dispersion is negligible in this limit and (26) reduces to 



dtH + dxH + 2Hd^H + ^dxxH + h 

6 5 



duH + ^dtxH + )^dxxH 



-dxxxxH = . (28) 



By making use of the first-order equivalence of the time- and space-derivatives of H, dtH 
—dxH, (28) reduces to the Kuramoto-Sivashinsky (KS) equation: 



dtH + dxH + 2HdxH + 



-6 + (3 



dxxH+ ^dxxxxH -0. 



(29) 



To look for the TW solutions to (29) in their moving frame, ^ = x — ct, we rescale the 
velocity as c = 1 + C, and the amplitude a.s H = eA and stretch the moving coordinate as 
C = BX: 

(30) 



-303"^ A + SeB'^A^ + B^T—A + —-A = 

dX dX'^ 



where T = |(5 + /3 and the condition lim^_^-|-oo H = has been enforced yielding a vanishing 
integration constant. Balancing each term in (30) gives CB^ ~ sB^ ~ B^T ~ 1 so that 
B ~ E ~ T3/2 and C ~ T^/s. 
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Writing B = T~^/^, e = fT^/^ and C = fiT^/yS then leads to 

-M + 2A^ + -A + ^A^0, (31) 

which is an ordinary-differential equation governing the TW solutions to the KS equation 
propagating at speed /i. Equation (31) admits a one-hump solitary-wave solution for a 
particular value of /i = /iq ~ 1.216 and an amplitude Amax ~ 0.784. Consequently, in the 
limit a — > 0, the speed c and the amplitude /imax of the one-hump solutions to (3) follow 
power laws of the form: 

c^l + 0.405 f^^^^ ^ I ^ 0.523 T^/^ . (32) 

Recall that the relations (32) have been obtained under the assumption of small amplitude, 
^max — 1^1, hence T must be small, which implies that both inertia and azimuthal 
curvature effects must be small. Viscous drag and gravity acceleration are the dominant 
physical effects. Following [9], we may refer to this regime as the drag- gravity regime. 

It is possible to interpret this regime as one characterized by the arrest of the instability 
growth by the flow advection. The dispersion relation corresponding to the KS equation 
(29) and governing infinitesimal perturbations around the base Nusselt flow of wavenumber 
k and angular frequency u is identical to (7) when /3 is substituted with T. One can then 
follow the same line of reasoning in going from (7) to (8) and define the ratio of the typical 
time of advection of the structure over its length, Tq = k^^, to the typical time of growth, 
Tg = [max(ci;j)]~^, of the instability: 

Therefore T ^ 1 corresponds to ^ Tg. The flow advection is much faster than the 
instability and the ratio Ta/rg controls the amplitude and speed of the observed structures 
as reflected by (32). 

2. Soliton-like regime 

We consider in this section the limit of a thick fllm, a ^ 1, though in practice, if ce 
is large, the axisymmetry of the flow might be difficult to achieve. We thus have (pia) ~ 
(3/4)alog(a), Ck ~ 4/35"^ and /3* ~ (3/4)^/3 (aGb)"^/^ 6 ~ fa log(a)r3/2 (aGb)"/^ while 
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in all cases 77 = (aGb)^/^. Therefore the conditions 5^1 and /3* ^ 1, necessary to sustain 
the weakly nonlinear analysis, can be justified if e.g. aGo = 0(1), i.e. rj = 0(1) and 
r3/2 For silicon oil vlOOO, T = 0.1 so that T^/^ ^ q.OS. [We note that strictly speaking 
the derivation of the WRIBL model (3) demands that F is at least 0(1) since the long- wave 
approximation is sustained only when surface tension effects are strong [17].] 
Neglecting inertial effects and looking for the dominant terms in (26) leads to 



6a 6a 6ari 6a 



3 1 
-^dxxtH + —dxxxH 
2 a 



. (34) 



-V log(") 

The first-order equivalence of the time and space derivatives of H, i.e. dfH = —{4:/{6a)dxH+ 
0(e£:^,e^£:), can again be utilized. Equation (34) is then simplified into 

dtH + ^d,H + ^Hd,H + -^d^^H + l^^$^d,,,H + ^d,,,,H = , (35) 
6a 6a 6arj a 6a 

the 'Kawahara' equation that was scrutinized numerically by Kawahara and coworkers [35, 
36]. The equation is also often referred to as the 'generalized KS' equation [20, 37, 38], and 
it is the KS equation appropriately extended to include dispersion (dxxxH). As before, we 
look for TW solutions in their moving frame, ^ = x — ct. Stretching the moving coordinate 
as = y/V^^ the amplitude as H = r]~^^'^A/2, and the speed as c = (4/3)a~-^(l + /i?7~^/^/4), 
then gives 

_,^ + 2^. + _^^ + fc_^^ + _^^.0, (36) 

with Sk = 3r/3/2 j^g^^ 

a) is a parameter that measures the relative importance of dispersion. 
Equation (36) governs TW solutions of the Kawahara equation propagating at speed /i. 
Since a ^ 1 and t] = 0(1), the dispersion parameter 6k is large, in which case the speed 
of the one-hump solitary waves solutions to (36) is given by /i ~ 0.32565i^. Solitary wave 
solutions of (3) in the limit a ^ 1 and 6 <^ 1 should therefore satisfy: 

— ^ 1 + 0.241og(a). (37) 

The speed of the one-hump solitary wave solutions to (3) is compared to the asymptotic 
limit (37) in Fig. 9(a). However, our computations show that for moderate values of a the 
wave speed is in fact closer to c/c^ ~ 1.65 + 0.24 log(a). The reason for this discrepancy can 
be traced in the violation of the assumption of small amplitude deviations from the base 
state (/i — 1 ^ 1) that was necessary to obtain (35) from (3). Indeed, TW solutions to the 
Kawahara equation have a speed and an amplitude that diverge as 6k becomes large. 
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(a) (b) 

FIG. 9. (a) Speed of solitary wave solutions to (3) when inertia is neglected {5 — )• 0) as function of 
the aspect ratio a. Dotted lines refer to (37) and to c/cfc ~ 1.65 + 0.241og(d). (b) Profile of the 
solitary wave solutions to (3) with (solid line) and without (dashed line) inertial terms (castor oil 
r = 0.45, R = 0.31 mm, 6 = 0.5, r] = 0.45 and a = 4.06). 

D. Speed to amplitude relation 

The weakly nonlinear analysis presented in subsection IV C provides us with the depen- 
dencies (32) of the speed and maximum height as functions of T = |5 + /3 for small film 
thicknesses. Combining these two relations give a linear dependence of the deviation ampli- 
tude /imax — 1 to the deviation speed c/ck — 1, where the wave speed c has been normalised 
with the kinematic wave speed c^: 

hra.. - 1 ^ 1.29 (^^ - 1^ (38) 

A Unear dependence of the speed as a function of amplitude was initially found by Chang [39] 

^max - 1 ~ - - 1 (39) 

by utilizing a normal form analysis of the TW solutions of the KS equation (29). This linear 
dependence is a characteristic of the drag-gravity regime and must be contrasted with the 
experimental relation 

/i^ax- 1-1.67(^^-1^ (40) 

obtained by Tihon et al. [40] for solitary waves running down a plane inclined at an angle 
5°. A linear dependence of /imax with respect to c/ck has also been found experimentally in 
the recent study by [5] for solitary waves on a relatively thick fibre. 
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FIG. 10. Deviation wave amplitude /imax — 1 versus relative wave speed c/c^ — 1 for the solitary- 
wave solutions to (3) for Go = 1.01, 0.236 and 0.168. The dashed lines correspond to the estimates 
(38) and (41). 

Similarly, in the drop-like regime, a law relating c and /imax can easily be obtained by 
recognizing that the wave amplitude /imax — 1 should be proportional to the distance hn — 1 
separating the fixed points. The constant of proportionality is determined by continuity 
with (38) in the limit /in — 1 <^ 1 for which hu — 1 ^ c/cfe — 1 (see Fig. 1). Next, by 
approximating hu with — 1/2 + A/3|c7c^^^^^nr74j we arrive at: 



Figure 10(a) presents the dependency of /imax — 1 on the relative speed c/cfc — 1 for the results 
in figure 2 corresponding to the WRIBL model (3) and to silicon oil v50. The direction of 
increasing reduced Reynolds number 5 along the curves is indicated by arrows. There is 
good agreement with relation (38) in the limit c ^ c^. A good agreement with relation 
(41) is also found in the case of the fast waves (c/cfc > 6), typical of the drop-like regime 
observed at small Goucher numbers. In the drag-inertia regime, however, at large values of 
5, the speed of the waves saturates to the limit Ccrit(<5, P/^) and no univocal mapping of the 
amplitude to the wave speed is observed. 

When the drop-like regime is affected by inertia for weakly viscous liquids, the linear 
relation (41) is no longer applicable reflecting the influence of the K mode on the drop-like 
waves and in particular the steepening of the wave fronts (see Fig. 10(b)). 

It is then clear from the above, that equation (41) provides a good approximation to the 




(41) 
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speed-to-amplitude dependence for solitary waves both in the drop-like and the drag-gravity 
regime. 

V. TRAVELING WAVES 

In this section, we compare periodic TW solutions of the WRIBL model (3) — i.e. limit 
cycles of the dynamical system (13) — with (77 7^ 0) and without viscous dispersion (77 — > 0) 
and of the CM equation (5) (both 6^0 and 77 —t- 0), for the experimental conditions 
considered by [2] and [4]. We revisit the description of the TWs branches of solutions 
given in [17] with reference to the characteristics of TWs of infinite extension, i.e. as TWs 
approach homoclinicity. In the experimental set-ups of [2, 4] the flow rate was controlled and 
maintained at a constant value at the inlet. Amplification of the ambient noise resulted in 
a wavy dynamics with waves traveling with constant shapes and speeds over long distances. 
Periodic TWs can be produced experimentally by applying periodic perturbations at the 
inlet [5, 41]. If the signal remains periodic in time at each location in space, an integration 
in time of the mass balance (3a) shows that the time average of the flow rate, qdt, 
where T is the period, is conserved all along the fibre and is equal to its value at the 
inlet which gives the condition qdt = 1/3 [42]. Periodic TW solutions of the model 

equations aiming to describe the experimental conditions must therefore verify the condition 
(q) = k/{2n) j^"'''qdi = 1/3, where k again denotes the wavenumber. 

Let us first consider the experimental conditions corresponding to the very viscous fluid 
considered by [2] (castor oil, F = 0.45; see Table I). The bifurcation diagram of TW solutions 
of the WRIBL model (3) and of the CM equation (5) are compared in Figure 11(a). The 
parameters in the figure correspond to regime 'c' reported by Kliakhandler et al. (gN = 
21 mg/s and R = 0.25 mm). We have normalised the wavenumber k with the reference 
knp = v^/ (1 + '5) corresponding to the marginal stability condition (real k ad u) for linear 
waves solutions to the CM equation k-^p corresponds to a dimensional wavelength 27i{R+hf^) 
proportional to the diameter of the liquid cylinder. 

Figure 11(c) shows corresponding wave profiles with regularly spaced streamlines in the 
moving frame, i.e. iso-contours of the function ijj = J^Ux rdr + c{R'^ — r'^) /2 at levels nqo/N, 
1 < n < N with = 10. The surface of the fibre corresponds to ip = and the free 
surface to ip = go- In the case of the WRIBL model (3), only one branch of TW solutions 

31 




(c) 

FIG. 11. TW branch of solutions bifurcating from the marginal stability conditions. (a) Normalised 
speed c/cfc as function of the normalised wavenumber k/k^ip where /crp = + a) (see text). 

Solid and dashed lines refer to (3) and to the CM equation (5), respectively, (b) Maximum thickness 
/imax (dashed line) and substrate thickness hs (dashed-dotted line) as functions of the normalised 
wavenumber k/k-^-p. The solid line refers to the relative maximum thickness h^i^^/hg. Labels 1 refer 
to model (3), whereas labels 2 refer to the CM equation (5). (c) Wave profiles and streamlines 
in the moving frame for solutions indicated by crosses in panel a; left: solutions to (3); right: 
solutions to (5). Parameters correspond to the experimental conditions of [2]: = 5.3 mg/s and 
R = 0.25 mm {Go = 0.138, a = 2.04, 5 = 0.01, ry = 0.18 and 13* = 2.1). 

has been found emerging from the marginal linear stability conditions at /c ~ /crp through 
a Hopf bifurcation, whereas a secondary branch has also been found by period doubling 
for the CM equation (5). We note that weakly nonlinear TW solutions of (3) travel at 
a smaller speed than the TW solutions of the CM equation (5) since the speed of linear 
kinematic waves is significantly affected by the second-order viscous effects [17]. At small 
wavenumbers, TW solutions of the branch emerging from k ~ /crp accelerate, become more 
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and more localized and terminate into the single-hump solitary waves discussed in § III. The 
speed, amplitude and shape of the solutions of (3) and (5) are comparable in this limit. 
Following the terminology introduced by Chang and co-workers for falling films on a planar 
substrate [43, 44], we refer to this branch of solutions as the '72 waves'. The branch of 
solutions to (5) emerging through period doubling terminates by slow waves with a shape 
made of a trough followed by capillary ripples. Following again [43, 44] for planar substrates, 
we refer to these TWs as the '71 waves'. 

For the conditions of Fig. 11, the fibre radius is small compared to the capillary length 
{Go = 0.138) and the Nusselt film thickness is comparable to the fibre radius (a = 2.04). As a 
consequence, the RP instability mechanism is strong and TWs have amplitudes comparable 
to the fibre radius, which corresponds to the typical situation of the 'drop- like' regime of the 
solitary waves discussed in § IV A. Yet, a direct transposition of the results obtained in § IV A 
to solitary-like wavetrains in the limit A; — is not possible since the reference thickness 
of the substrates of the solitary-like waves is not constant (recall that in the treatment 
of solitary wave solutions, the constant film thickness far from the solitary waves was the 
reference thickness). We have computed the parameter /3* based on the substrate thickness 
hg defined by the position of the fixed point in the phase space when the corresponding 
limit cycle approaches homoclinicity. The maximum of the relative thickness /imax/^s, 6.5 
here, is reached for a value of the local number /3* close to its maximum, namely 3.65, as 
expected since the waves' characteristics are piloted by the balance between the advection 
by the flow and the RP instability in the drop-like regime. Figure 11(b) shows the evolution 
of the maximum thickness, /imax? of the substrate thickness hs and of the ratio of the two 
as function of the normalised wavenumber k/kup. At small k, a TW is solitary-like and 
travels on a portion of nearly flat substrate of increasing length. The mass transported by 
the wave decreases in comparison to the mass carried by the substrate and hs asymptotes to 
the Nusselt film thickness. As a consequence hs increases as k tends to zero. This trend is 
followed by the maximum height /imax- As a TW gets more and more localized, it tends to 
have bigger and bigger size. However, the maximum relative thickness, /imax/^s? is a rapidly 
decreasing function of k at small wavenumbers. Indeed, larger substrate thicknesses imply 
weaker RP instability mechanisms and stronger advection: /3* (not shown) follows a trend 
similar to /imax/^s- We therefore obtain a rather paradoxical picture: Despite the weakening 
of the RP instability — thus the lowering of /imax/^s — , a TW has its absolute amplitude 
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augmented as k tends to zero. 

For the experiments by Kliakhandler et ai, the RP instabihty mechanism is strong while 
inertia and viscous dispersion are weak (5 = 0.01, 77 = 0.18) which explains the good agree- 
ment obtained with the results from the CM equation (5). However, as the RP instability 
is weakened by raising the Goucher number, streamwise viscous dispersion should play an 
increasingly important role. We have checked this by computing the TW solutions of (3) 
and (5) corresponding to a flow rate q-^ = 10 mg/s and a larger radius R = 0.5 mm, i.e. for 
6 = 0.01, f3* = 1.16 and rj = 0.22, parameters that are comparable to those corresponding to 
the regime 'c' discussed by Kliakhandler et al. {5 = 0.01, /3* = 2,1 and 77 = 0.18) but with 
a value of (3* ~ two times smaller. Branches of solutions are compared in figure 12(a) in 
the (/c/Zcrp, c/cfc)-plane. 72 fast TW solutions of (5) are again observed to emerge from the 
marginal conditions k ~ /crp whereas 71 solutions are obtained through a period doubling 
bifurcation of the 72 branch. We note again that only the 72 branch of solutions can be 
found for the WRIBL model (3), with no period doubling bifurcation being detected in this 
case. We therefore conclude that the inclusion of second-order viscous effects reduces the 
number of wave families and, as a consequence, it may drastically simplify the complex se- 
quence of bifurcations and topological structure of associated bifurcation diagrams. In fact, 
such a reduction of the number of wave families by the second-order viscous terms was also 
evidenced in the planar case [45]. 

In comparison to Fig. 11(a), the 72 branches portrayed in figure 11(a) now deviate signif- 
icantly from each other even at small wavenumbers. The TW solutions of model (3) have a 
larger amplitude than the TW solutions of equation (5) as indicated in Fig. 12 which com- 
pares the wave profiles of waves of equal wavenumber. When approaching homoclinicity, 
this effect is enhanced, and solutions of (3) travel with a notably higher speed and larger 
amplitude than solutions of (5). They are also preceded by smaller capillary ripples. 

Figure 12(b) compares hg and /imax- Again, the substrate thickness and the maximum 
thickness have the same trend as k tend to zero whereas the maximum relative thickness 
hraax/hs cvolvcs in the opposite direction. However, a comparison of Fig. 12(b) to Fig. 11(b) 
reveals that the maximum reached by /imax/^s is ~ three times smaller than with the ex- 
perimental conditions considered by [2] . Computations of the local values of the parameters 
based on the substrate thickness show that f3* does not exceed 1.41 and decreases with k 
whereas rjs increases. The local parameters are /3* = 1.2 and rjs = 0.19 for the wave profile 
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FIG. 12. TW branch of solutions bifurcating from the marginal stabihty conditions. Parameters 
correspond to castor oil (F = 0.45), = 10 mg/s and R = 0.5 mm {5 = 0.01, rj = 0.22, a = 1.15 
and /5* = 1.16). See caption of figure 11. 

labelled (iii) in Fig. 12 so that viscous dispersion effects start to be comparable with the 
RP instability mechanism. We conclude that viscous dispersion has a non-trivial amplifying 
effect on the RP instability mode. 

For less viscous fluids like Rhodorsil silicon oil v50, one may expect a more significant 
effect of inertia and of the K mode of instability. In Fig. 13(a), we compare the substrate, 
maximum and relative maximum thicknesses, hg-, /^max and /imax/^s, respectively, corre- 
sponding to the TW branch of solutions to (3). Parameters are chosen to correspond to 
one of the experimental conditions considered by [4] [R = 0.23 mm, = 151 mg/s) and 
to a Goucher number Go = 0.155 close to the corresponding one for the experiments in [2] 
(Gb = 0.137). 

As k is lowered and the TWs approach homoclinicity, hs and /imax grow, a trend similar 
to what is observed in figures 11(b) and 12(b) for a more viscous fluid. However, the relative 
amplitude /imax/^s follows a surprising non-monotonic behaviour with an S-shape curve in 

35 



(a) (b) 

FIG. 13. (a) Maximum thickness /imax (dashed line) and substrate thickness hg (dashed-dotted 
line) versus the normalised wavenumber k/kup. The solid line refers to the relative thickness 
/iinax/^s- (b) Local reduced Reynolds number 6s and /3* versus k/k^p. Parameters correspond 
to Rhodorsil v50 silicon oil (F = 5.48), = 151 mg/s and R = 0.23 mm {Go = 0.155, 6 = 3.9, 
a = 3.0 and r] = 0.36). 

the plane [h^k/k-^-p)- This behaviour can be understood by examining the local values of 5 
and /3* based on the substrate thickness and presented in figure 13(b). As k is lowered 
to zero, /3* decreases whereas 5s varies in the opposite direction. As a result, the two curves 
ultimately cross and very long waves correspond to a relatively stronger K instability mode 
than the initially dominant RP instability. We thus recover the usual trend observed with 
the solitary waves propagating on a planar film: Larger substrates imply larger amplitudes 
/imax and also larger relative amplitude /imax/^s [44, 46]. 

We now turn to the experimental conditions considered in [5] corresponding again to 
Rhodorsil silicon oil v50 but with a larger fibre, R = 0.475 mm. The chosen flow rate 
gN = 296 mg/s gives values of the reduced Reynolds number S and of the viscous dispersion 
parameter t] which do not vary significantly (compare e.g. 6 = 4.1 and t] = 0.44 to 6 = 3.9 
and 7] = 0.36). Inertia and viscous dispersion effects are therefore comparable whereas the 
RP instability mechanism is significantly lowered (/3* is reduced from 1.30 to 0.77). As for 
R = 0.5 mm and = 10 mg/s, we found a single branch of TW solutions bifurcating from 
the marginal stability curve and the obtained bifurcation diagram (not shown) is similar to 
Fig. 12. Figure 14(a) compares the corresponding substrate thickness hs, the absolute and 
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FIG. 14. (a) Maximum thickness /imax (dashed hne) and substrate thickness hg (dashed-dotted hne) 
versus the normahsed wavenumber fc/ZcRp. The sohd hne refers to the maximum relative thickness 
hmax/hg- (b) Local reduced Reynolds number 6s and /3* versus k/k-^p. Parameters correspond to 
Rhodorsil v50 silicon oil, (/n = 296mg/s and R = 0.475 mm (Go = 0.32, 5 = 4.1, a = 1.7 and 
r] = 0.44). 

relative amplitudes /imax and h^i^^/hs- As k is lowered, all curves exhibit a similar monotonic 
shape similar to what is expected in the planar case. Again, this can be understood by 
examining the local values of the two governing parameters (cf. Fig. 14(b)). The local 
reduced Reynolds number /3* does not exceed unity whereas 6s sharply increases for very 
long waves, which signals the prevalence of the K mode over the RP mode. 

Interestingly, TWs can be found at a wavenumber k higher than the critical wavenumber 
kc ~ ^Rp corresponding to the marginal stability conditions, in which case the Nusselt 
solution is linearly stable. This subcritical bifurcation of TW solutions from the Nusselt 
uniform film occurs at small radii of the fibre, i.e. small Goucher number (see figures 11 
and 12). The subcritical onset of TWs has also been reported in the recent work by Novbari 
and Oron [47] based on an 'energy integral method' though these authors did not give a 
physical interpretation of this phenomenon. Since the RP instability mechanism is strong 
when Gb ^ 1, we can conjecture that the onset of subcriticality is related to surface tension 
effects associated with the azimuthal curvature (hence this effect is absent in the planar 
case). From a thermodynamic viewpoint, surface forces tend to reduce the free energy of 
the system and equilibrium is reached when the contact area between the liquid and the 
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FIG. 15. Traveling wave solutions to the CM equation (5) for a = 8 and a = 0.5. Thick and 
thin solid lines refer to Gb = 0.02 and Go = 0.055, respectively. Curves labeled 1 (respectively 2) 
correspond to ci = 8 (q = 0.5). 

surrounding gas is minimum. Obviously, the flow is out of equilibrium and a thermodynamic 
argument must be taken with care [48]. Yet, when 5 ^ 1, 77 ^ 1 and Gb ^ 1, model (3) 
reduces at leading order to the long-wave Young-Laplace equation and the shape of the TWs 
thus only slightly differs from the shape of static drops on fibres as was shown in § IV A and 
in Fig. 6. 

We have computed the TW solutions to the CM equation (5) for Go = 0.02 and Go = 
0.055, and for a = 0.5 and a = 8. The maximum amplitude /imax, and the interfacial area 
A of the waves normalised by the area of the uniform film solution of the same volume 
are displayed in Fig. 15. The subcritical onset of the TW solutions is accompanied by a 
normalized interfacial area lower than unity, and therefore a lowering of the free surface 
energy in comparison to the Nusselt uniform film solutions, which supports our conjecture 
that subcriticality is promoted by capillary effects. We note that the minimum wavenumber 
at which TWs are observed seem to depend only on a independently of Go. Conversely, 
at a given value of a, the amplitude of the waves is strongly affected by the value of Go, 
a small Goucher number implying larger waves as expected, since the saturation number 
l3 = a'^^^Go~^^^ is also larger. 

In Fig. 16 the maximum of the normalised wavenumber k/kup of the TW solutions of 
the CM equation is depicted as function of the aspect ratio a for Go = 0.02 and Go = 
0.055. We note again that the location of this local maximum only weakly depends on the 

38 



Goucher number. The shape of the corresponding TW solutions (not shown) is also nearly 
symmetrical. We thus conclude that the subcritical behavior of the TW solutions is not 
affected by gravity effects at small values of the Goucher number but is strongly dependent 
on the parameter a and therefore on the geometry of the base flow. 

By neglecting gravity effects and assuming a wettable substrate, the shape of a static 
drop is governed by its length and volume. Therefore, for a normalized area A set to unity, 
the shape of a static drop depends only on the aspect ratio a of the uniform film solution 
with the same volume. The solid line curve in figure 16 represents the wavenumber of static 
drops with unit normalized areas. Above this curve, A> \ and the static drop solution is 
not energetically favorable, whereas below it the static drop solution is favored. Indeed, large 
drops have a nearly spherical shape with an interfacial area that is lower than the uniform 
film with the same volume, with further lowering of the interfacial area being forbidden 
by the presence of the fibre. For a given amount of liquid, the thicker the uniform film 
solution is, the lower the normalized area of the static drop would be. Therefore, the range 
of energetically favored wavenumbers of the static drop solution increases with the aspect 
ratio a as can be observed from figure 16. 

We can now deduce two more arguments from Fig. 16 which support a capillary origin 
of the subcritical behavior of the branch of TW solutions that emerge from the Hopf bifur- 
cation of the uniform film solution: (i) Since all points lie below the solid curve, TWs only 
exist when an energetically favorable static drop solution is available; (ii) The trend of the 
maximum of the TW wavenumber with respect to the aspect ratio a is similar to the trend 
of the boundary separating energetically favored and unfavored static drop solutions. 

VI. PHASE DIAGRAM 

We are now in a position to give a phase diagram of the different regimes found for all 
possible TW solutions on a film flow down a fibre. The onset of the 'drop-like' regime and 
the 'drag- inertia' regime corresponds roughly to /3* ~ 1 and 5 ~ 1. The 'soliton-like' regime 
arises when the instability mechanisms are weak (5^1 and /3* ^ 1), the film is thick, 
a = 0(1), and viscous dispersion is strong, t] = 0(1). Finally, the 'drag- gravity' regime 
is observed when all other effects are weak (5^1 and /3* ^ 1 and ?7 <^ 1). Therefore, 
a phase diagram can be obtained for a given fluid, thus a given Kapitza number F, by 
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FIG. 16. Location of the maximum of the normahsed wavenumber k/k-^p as function of the aspect 
ratio Q for the TW solutions of the CM equation (5). Crosses (respectively squares) correspond to 
Go = 0.02 {Go = 0.055). The solid line denotes the static drop solution to (Bl) with a normalized 
area A = 1 (see text and Appendix B). 

drawing the curves S = 1, r] = 1 and /3* = 1 in the plane (a, Go). Since r] = (aGb)^/^ 
and /3* = {ack{a) /[Go'^{l + a)^]}^/^ are functions of a and the Goucher number only, the 
corresponding curves /3* = 1 and t] = 1 are independent of the working fluid considered. 
Thus, 5 = 1 is the only boundary that moves in the plane [a, Go) when F is varied. Figure 17 
is a tentative representation of the phase diagrams for the four working fluids considered in 
this study, from weakly viscous fluids like water with a high Kapitza number, F = 3376, to 
highly viscous fluids like silicon oil vlOOO corresponding to a small Kapitza number, F = 0.10. 
In the case of water, a large overlap region exists for the 'drop-like' and the 'drag- inertia' 
regime corresponding to a mutual reinforcement of the two K and RP instabilities. The 
'soliton-like' regime takes over at relatively high viscosities where the curve 6 = 1 moves 
below the curve rj = 1. 

VII. CONCLUSIONS AND DISCUSSION 

We have presented new results and insights on the characteristics of axisymmetric waves 
propagating down a fibre. Our analysis was based on the two-equation model derived in 
[17] using a weighted-residuals approach, equation (3). The model accounts for all physical 
effects, namely inertia, azimuthal curvature and viscous dispersion and has been validated 
in the studies by Ruyer-Quil et al. [17] and Duprat et al. [5] through direct comparisons to 

40 




1 1.2 



(a)water (r = 3376) 




(c)castor oil {T = 0.45) 




1 1.2 



(b)silicon oil v50 (F = 5.48) 



4 
3.5 

3 
2.5 
ifi 2 
1.5 
1 

0.5 






y? =\ DI 
























\/3* = 1 






) DG 

, — III 





0.2 0.4 0.6 O.i 

Go 



1 1.2 



(d)silicon oil vlOOO (T = 0.10) 



FIG. 17. Maps of the different regimes in the Go — a parameter space for fluids of increasing 
viscosity. The different curves are the loci of 6 = 1, r] = 1 and /?* = 1. 'DI' refers to the 'drag- 
inertia' regime, 'DL' to the Rayleigh-Plateau regime, 'SOL' to the 'soliton-like' regime and 'DG' 
to the 'drag-gravity' regime. 

the experiments by Kliakhandler et al. [2] and Duprat et al. [4, 5]. 

We first focused on isolated waves running on a constant thickness substrate, or solitary 
waves. The dynamics of the film seems to be dominated by these structures even when the 
flow becomes disordered [2, 3]. We examined in detail, via both asymptotic analysis and 
through elements from dynamical systems theory, the shape, speed and amplitude of the 
waves for four fluids of increasing viscosities: Water, Rhodorsil silicon oils v50, vlOOO and 
the castor oil utilised in the experiments by [2]. 

We identified four distinct regimes corresponding to the competition of the two instabil- 
ity modes, the K and RP modes, prompted by inertia and azimuthal curvature (the fibre 
curvature effectively), respectively, with the viscous dispersion (i.e. the second-order axial 
viscous diffusion and second-order viscous contributions to the tangential stress at the free 
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surface) and with the advection of the structures by the flow, which results from the balance 
between gravity and viscous drag. Two of these regimes are similar to what is found in the 
planar geometry [9]. The 'drag-gravity' regime corresponds to the predominance of the flow 
advection over the instability mechanisms, either when inertia effects are weak, i.e. for 5 ^ 1 
or when the azimuthal curvature effects are non-dominant, /3* ^ 1. In both cases it is possi- 
ble to interpret the 'drag-gravity' regime as one where the instability growth is arrested by 
the flow which determines the amplitude and speed of the solitary waves as reflected by the 
asymptotic relations (32). The 'drag-inertia' regime is observed at large reduced Reynolds 
numbers, 6^1, when the wave characteristics are determined by the balance of inertia, 
drag and gravity. We have obtained the asymptotic limit of the speed and showed that the 
rate of convergence to this limit is governed hj rj/ 6"^ cc Ee~^. 

The 'drop-like' regime corresponding to the predominance of the RP instability mecha- 
nism over the flow advection. It is specific to the cylindrical geometry and is observed for 
small fibre radii R compared to the capillary length Z^, that is at small Goucher numbers 
Gb, when the typical time of growth of the RP instability is greater than the advection time 
of a wavy structure, i.e. for /3* ^ 1. The maximum reachable amplitude and speed of the 
waves in this regime is governed by the radius R of the fibre and the balance of gravity and 
viscous drag. Comparisons to quasi-static drop solutions of the Laplace- Young equation 
(Bl) sliding down a fibre with a speed verifying the Landau-Levich-Derjaguin law (B4) 
show excellent agreement, even in the case of spherical drops where the long- wave lubri- 
cation assumption does not strictly apply. In this regime, waves have a drop-like nearly 
symmetrical shape determined by capillary effects. The thickness of the substrate film on 
which the drops slide is governed by the balance of viscosity and capillarity. 'Drop-like' TW 
solution branches subcritically emerge from the Nusselt uniform film branch. We have given 
an explanation for this subcritical onset based on geometric and thermodynamic arguments 
and thus completed its recent investigation in [47]. This phenomenon arises from capillary 
effects and depends only on the aspect ratio a for sufficiently low Goucher numbers. 

We have also found a possible fourth regime for very viscous fluids and thick films a = 
0(1)), for which both K and RP instability mechanisms are weak (5 ~ 1 and /3* ^ 1) and 
viscous dispersion is significant (77 = 0(1)). This 'soliton-like' regime corresponds to the 
balance of the nonlinearities with the dispersion induced by second-order viscous effects, 
with the speed and amplitude of the solitary waves being functions of the logarithm of the 
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aspect ratio a. 

Our study of the solitary-wave solutions has been followed by construction of the TW 
branches of solutions corresponding to the experimental conditions for which the average 
flow rate is the true control parameter, with the substrate thickness being determined by 
the solution itself. If the substrate thickness and the maximum amplitude /imax grow 
when TWs approach homoclinicity, the ratio of the two, h^iay^/hg, evolves in a manner that 
strongly depends on which instability mechanism is dominant. If the RP instability is 
dominant, /imax/^s decreases as the wavenumber k tends to zero, whereas if the K mode is 
dominant h^g^^/hs has an opposite trend. This picture can be even more complex since the 
predominance of the instability modes can be exchanged by varying the periodicity of the 
waves and thus the substrate thickness (cf. Fig. 13). The selected wave regime depends not 
only on the properties of the Nusselt flow at the inlet but also on the periodicity selected by 
the system. Indeed, the boundaries separating the different regimes are not only functions 
of Go and a but also functions of the thickness of the substrate, which is determined by the 
typical distance separating solitary-like waves. Therefore, the phase diagrams displayed in 
Fig. 17 must be taken with caution. The wave selection process of a noise-driven falling film 
is the complex result of the linear amplification of inlet perturbations and the downstream 
nonlinear interaction mechanisms [21, 49, 50]. 

Noteworthy is that in our previous study on the fibre problem [17], we have shown that 
the wave selection process of the noise-driven wave dynamics down the fibre is strongly 
affected by viscous dispersion. This effect is generally weak in the experiments devoted 
to falling films on planar substrate, where working fluids were often weakly viscous, i.e. 
water [40, 41, 51]. In the case of films down fibres, the oils used in the experiments are much 
more viscous than water which explains that viscous dispersion can be dominant and can 
promote a 'soliton-like' regime. In this regime and despite the dissipative nature of the flow, 
it is possible to observe the formation of solitons, i.e. solitary waves whose shape and speed 
are not altered by collisions with other solitary waves (such waves are still dissipative but 
they share several common features with solitons in conservative systems [52]). We explore 
such effects in Fig. 18 which depicts a simulation of the interaction between two solitons. 
Chosen parameters correspond to Rhodorsil silicon oil vlOOO and Go = 0.6. Two pulse-like 
perturbations of different amplitudes are initially placed in the computational domain (cf. 
panel a). The excess mass carried by these two perturbations are then drained behind the 
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FIG. 18. Simulation of the response of the film in the 'soliton-like' regime. Initial conditions 
consist of two pulse-like perturbations of different amplitudes. Parameters are 5 = 0.18, a = 1.9 
and r/ = 1.2 (Rhodorsil silicon oil vlOOO, = 300 mg/s and R = 0.89 mm), (a-d) Snapshots of the 
film thickness at increasing times; (e) Spatio-temporal diagram. Vertical and horizontal ranges are 
1.25 m and 28 s respectively. Elevations (depressions) of the free surface are coded in dark (light) 
grey. 



two pulses creating two wave-packets. The spatial extend of these packets grows in time 
due to the convective instability of the film [4, 17]. Panel b shows the interaction of the 
second pulse with the mass ejected by the first pulse. The spatio-temporal diagram shown 
in panel e indicates that the second pulse first overtakes this excess mass and next the first 
pulse without absorbing them (cf. panels c and d). The coalescence of the two pulses is 
accompanied by a phase shift but no notable modifications of the speeds and amplitudes 
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of the pulses as in Hamiltonian systems: The initial wave profiles get superimposed as the 
waves collide and reappear as the waves move apart. Yet, a slow evolution of the amplitude 
and speed of the two solitons can be observed towards the speed c = 7.5 cm/s and amplitude 
^max = 3.3 mm of the infinite-domain solitary-wave solution for the given set of parameters. 

Our hope is that this new evidence of existence of solitons on a liquid film flowing down a 
fibre, a dissipative system as opposed to a Hamiltonian one, may motivate a renewed interest 
for the experimental investigation of the resulting wave dynamics and, in particular, on the 
role of viscous dispersion. 
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Appendix A: Coefficients of the model (3) 



The expressions of the coefficients entering the averaged momentum balance (3b) consist 
of ratios of polynomials in b and log(6), where b = 1 + ah is the ratio of the total radius 
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R + h to the fibre radius R: 



4>=[3 ((4 log(6) - 3)6^ + 46^ _ i)] / [i6(6 - l)^] , (Ala) 
F = 3Fj[16{b-lf<l)F,], (Alb) 
Fa = -3016^ + 6226^ - 4416^ + 4 hg{b) {l97b^ - 2346^ + 6 hg{b) 

X [16 log{b)b^ - 36b^ + 226^ + 3] 6^ + 78^2 + 4} 6^ + 13052 _ 10 , (Ale) 
Ff, = 17b^ + 12 \og{b) [2 log{b)b^ - 36^ + 2] 6^ - 306"^ + 156^ - 2 , (Aid) 
G' = ay[64(6-l)Vi^b], (Ale) 
= % |4 log(6) [-2206^ + 4566^ - 3036^ + 6 \og{b) {61b^ - 6%^ 

+4 Iog(6) (4 log{b)b^ - I2b* + 76^ + 2) 6^ + 952 + 9) 52 ^ 58^2 + 9] 52 
+ (6^ - 1) ' (1536^ - 1456^ + 536" - O } ' (Alf ) 

/ = 64(6-l)50V[3i^6], (Alg) 
J = 3J,/[128(6-l)Vi",], (Alh) 
Ja = 9 I (4906^ - 2056^ - 2356^ + 736^ - 3) (6^ - l)^ 

+46^ log(6) [2b^ \og{b) (72 log(6) (2 log(6)6^ - 66^ + 6^ + 6) 6^ 
+ {b -l){b+ 1) (5336^ - 1096^ - 4516^ + 15)) 
-3 (6^ - 1)^ (1876^ - 436^ - 1346^^ + 176^ + l) } , 
K = 3Kj[m\b-lf(PF^], 

Ka = 4bHog{b) (2336^ - 3606^ + 121og(6) (l21og(6)6^ - 256^ + 126^ + 9) b^ 
+546^ + 88fe2 - 15) - (^2 - 1)2 (211b^ - 1346^ - 566^ + 306^ - 3) , 
L = La/[8b{b-iy<j)Ft], 

La = 4b^ log(6) {6 log(6) (12 log(6)6^ - 236^ + 186^ + 3) + (b - l){b + 1) 

X (956^ - 796^ - 76^ + 3) } - {b^ - if {82b^ - 77b^ + Ab^ + 3) , 
M = 3 + [24 log(6)6^ - 256^ + 486^ - 366^ + 166^ - 3] /[2b'^Fb] . 



(Ali) 
(Alj) 

(Alk) 
(All) 

(Aim) 
(Aln) 



We note that in Appendix B in [17] a small misprint can be found in the definition of the 
factor J — a missing factor of three — that is here corrected. 
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Appendix B: Static drops on coated fibres: Computations and asymptotic analysis 



The shape of an axisymmetric drop sitting on a vertical fibre has been computed numeri- 
cally by Kumar and Hartland [53] and determined analytically by Carroll [54] by neglecting 
gravity effects. When the contact angle of the liquid with the solid fibre vanishes, the ana- 
lytical solution corresponds to an unduloid [7, 55] that can be written parametrically using 
elliptic integrals of the first and second kind. However, the analytical solution is cumber- 
some to use and requires numerical evaluation of the different integrals involved. Hence, 
we choose to determine the solution by solving numerically the Laplace- Young equation 
parametrically rewritten as [53]: 

1 sin(/) 

(Bla) 



ds \ds J ^ Rt r 

— = cos0, and — = sin0, (Bib) 
ds ds 

where the length scale is the fibre radius, R. The radial and axial coordinates are denoted by 

f and X, respectively, s denotes the curvilinear arc length along the drop interface whereas 

(f) is the angular inclination of the drop interface to the radial axis whereas {d(f)/ds)^ + 1/Rt 

denotes the mean curvature at the top of the drop, s = 0. The set of equations is completed 

by the boundary conditions: 

= 7r/2, r = Rt, x = 0, at 5 = 0. (Blc) 

The contact area A separating the liquid and gas phases and the volume V of the drop has 
been computed by solving: 

^ = 27rf and ^ = nir' - R',) (Bid) 
ds ds 

We have solved system (Bl) using the Auto07p software [30]. We have adjusted Rt to the 
coated fibre radius 1 + a, and the drop volume V to the volume of the corresponding solitary 
wave, after subtracting the volume of the residual film. 

The speed of the quasi-static drops sliding on a vertical fibre can be estimated using 
the Landau-Levich-Derjaguin theory [56, 57]. For this purpose we divide a large amplitude 
sliding drop in two separate regions. The 'inner' one corresponding to the thin films at the 
upper and lower ends of the drop where viscosity balances capillary forces, and the 'outer' 
one, the quasi-static drop itself which is governed by the Laplace- Young equation given 
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above. In the inner region, h <^ R and the equation to be solved reduces to (6), whose TW 
solutions are governed by: 

— {l + (3h' + h"')-c{h-l)-l=0 (B2) 
3 

Following [13], we introduce the inner coordinates X = c~^^^^ where c ^ 1 and get to 
leading order 

h'" = 3{h - l)/h^ (B3) 

which is the so-called 'Bretherton equation' [13, 58], originating from Bretherton's work on 
the motion of a long gas bubble in a capillary tube. In the upper region, the solution of the 
Bretherton equation with boundary conditions /i = 1, /i' = /;," = at X = is monotonic. 
The numerical solution of (B3) yields limx^oo ^" = 1.34. The speed of the drops is then 
selected by asymptotically matching the solutions of (Bl) and (B3) in the upper overlap 
region between inner and outer domains. This can be easily done by imposing that the 
Laplace pressures corresponding to the two solutions are equal in this region giving: 

f).}^l. IM^. ,B4) 

as J ^ Kt a 

where Ca = CdmpsBo = nc^rops/c in the limit a <C 1. 
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